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A pair of quantum dots, coupled to each other through a point contact, can exhibit Coulomb 
blockade effects that reflect the presence of an oscillatory term in the dots' total energy whose value 
depends on whether the total number of electrons on the dots is even or odd. The effective energy 
associated with this even-odd alternation is reduced, relative to the bare Coulomb blockade energy 
Up for uncoupled dots, by a factor (1 — /) that decreases as the interdot coupling is increased. 
When the transmission coefficient for interdot electronic motion is independent of energy and is the 
same for all channels within the point contact (which are assumed uncoupled), the factor (1 — ,/") 
takes on a universal value determined solely by the number of channels A^ch and the dimensionless 
conductance g of each individual channel. When an individual channel is fully opened (the limit 
g — > 1), the factor (1 — /) goes to zero. 

When the interdot transmission coefficient varies over energy scales of the size of the bare Coulomb 
blockade energy Up, there are corrections to this universal behavior. Here we consider a model in 
which the point contact is described by a single orbital channel containing a parabolic barrier 
potential, with ujp being the harmonic oscillator frequency associated with the inverted parabolic 
well. We calculate the leading correction to the factor (1 — /) for N^h = 1 (spin-split) and Ach = 2 
(spin-degenerate) point contacts, in the limit where g is very close to 1 and the ratio 2TiUp/hujp is 
not much greater than 1. Calculating via a generalization of the bosonization technique previously 
applied in the case of a zero-thickness barrier, we find that for a given value of g, the value of (1 — /) 
is increased relative to its value for a zero-thickness barrier, but the absolute value of the increase 
is small in the region where our calculations apply. 

PACS: 73.23.Hk,73.63.Kv,71.10.Pm,72.10.-d 



I. INTRODUCTION 

In recent years, there have been a number of theo- 
retical and experimental studies of the manner in which 
Coulomb blockade effects in a metallic particle or semi- 
conductor quantum dot disappear when the conducting 
island becaaaes more electrically connected to its envi- 
ronment .Erllil Here we focus on a system in which two 
symmetric quantum dots are defined by applying nega- 
tive voltages to gate electrodes that lie on the surface of a 
semiconductor heterostructure above a two-dimensional 
electron gas (2DEG). We assume that the dots are joined 
by a quantum point contact containing a single orbital 
channel that is almost perfectly transmitting at the Fermi 
energy, but that the dots are isolated from their respec- 
tive leads by comparatively large tunnel barriers. In this 
geometry, information about the Coulomb blockade en- 
ergy Up involved in the transfer of electrons from one dot 
to the other can be obtained by observing the positions 
of Coulomb blockade peaks in the conductance across 
the entire system (from one lead to another) when that 
conductance is plotted as a fuiiction of voltages on gates 
coupled to each of the dots.ETEl 

Previous analyses of the disappearance of the Coulomb 
blockade for two-dot systemstd containing such a par- 
tially open point contact have characterized the contact 
by a number Nc\i of degenerate one-dimensional (ID) 
channels and by the dimensionless conductance g of each 



channel (where < 5 < 1).lI □ This characterization is 
complete in the limit where the electronic transmission 
amplitude through the contact is independent of energy 
for electron energies that differ from the Fermi energy 
by no more than an amount comparable to Up. The en- 
ergy independence of the transmission amplitude means 
that this characterization corresponds to a potential bar- 
rier which is sufficiently thin that it can be modeled as 
a delta function. Because the Coulomb blockade energy 
Up is much smaller than the Fermi energy Ep, working 
in this "delta-function barrier limit" can yield good re- 
sults even though the product of the barrier width and 
the Fermi wave vector kp \s generally much greater than 
1 (with kp being the value of the Fermi wave vector in 
the 2DEG far from the barrier). 

Nevertheless, it seems important to investigate further 
the consequences of relaxing the assumption of a delta- 
function barrier limit. For one thing, it is possible to 
generate wider barriers using an appropriate gate geom- 
etry, and one would like to understand at what point the 
delta-function-limit calculations break down. Secondly, 
estimating the corrections due to a barrier's finite thick- 
ness provides a valuable check on the delta-function- limit 
results. 

Another reason for interest is that recent experiments 
on transmission through quantum point contacts have 
shown unexpected structure (e.g., an apparent conduc- 
tance plateau near 0.7(2e^//i) at intermediate tempera- 
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turestS) whose origin is only poorly understood. Such 
results suggest that transmission through a point con- 
tact may have a nontrivial energy dependence, and such 
a dependence could well have an effect on the breakdown 
of the Coulomb blockade in the coupled-dot geometry. 
With the aim of explaining such effects in mind, it is a 
helpful first step to study theoretically the effects of en- 
ergy dependences in a simpler situation where the barrier 
potential is known and many-body effects are reduced to 
a bare minimum. In this vein we consider here a model in 
which the two dots are separated by a parabolic barrier 
of nonzero width and the electron-electron interaction is 
taken to be constant for any two electrons located any- 
where on the same quantum dot. 

This model presents challenges for "bosonization tech- 
niques" that characterize point-contact constrictions as 
one-dimensional fermionic seas whose low-energy degrees 
of freedom can be expressed in terms of bosonic den- 
sity and phase variables .Ej1l3 In such "bosonized" mod- 
els, the behavior of incompletely opened channels is 
commonly studied by introducing a zero-width delta- 
funfttinftjbarrier at a specific point in the one-dimensional 
sea.cHa'c3Eil By way of contrast, this paper seeks to show 
how a bosonized model of a one-dimensional system can 
describe the single-particle effects of replacing a delta- 
function barrier with a barrier that more realistically 
corresponds to a finite- length constriction — i.e., a barrier 
of nonzero width and therefore nontrivial single-particle 
transmission properties. 

Before proceeding, we should describe more fully what 
is aheady understood about systems of two symmetric 
quantum dots connected by a single orbital channel con- 
taining degenerate spin modes. In such systems, an en- 
ergy scale Up characterizes the energy cost of moving 
electronic charges between different segments of the sys- 
tem. Moreover, the system's Coulomb blockade behavior 
reflects the presence of an energy term proportional to 
Up that oscillates between one value held when the total 
number of electrons on the dots is even and ajjather value 
held when the total electron number is odd.ElB 

Previous work, both experimental and theoretical, has 
shown that as the conductance of the interdot point con- 
tact is increased, the energy scale associated with this 
even-odd alternation is reduced by a factor of (1 — /), 
where / goes to 1 when the channel is fully open. When 
the transmission coefficient for electronic motion within 
this channel is independent of energy (i.e., the delta- 
function barrier limit), the factor (1 — /) takes on a uni- 
versal value determined solely by the number of point 
contact channels iVch (the channels being assumed degen- 
erate and uncoupled) and the dimensionless conductance 
g of each individual channel (If G is the total conduc- 
tance of the point contact, g = G/iVch(e^//i).). We refer 
to the value / in the factor (1 — /) as the fractional peak 
splitting because it has been measured by observing the 
relative separation of conducta|iiGe|-,peaks in a series of 
Coulomb blockade experiments.u0u 

This paper studies the corrections to /, and therefore 



to (1 — /), that result when the assumption of a delta- 
function barrier limit is relaxed — in other words, when 
the interdot barrier is more realistically treated as hav- 
ing a finite height and nonzero width. A prior study 
of the effects of a nonzero-width barrier concentrated on 
the limit of weakly coupled dots (the limit g 0).E£l This 
study revealed that, if the one-dimensional channels be- 
tween two dots, or between a dot and a lead, contain 
a tunneling barrier of finite height Vq (in energy units) 
and of nonzero width ^, the behavior of such systems is 
responsive to another energy scale, W, that character- 
izes the energy range (for electrons incident on the bar- 
rier) over which the probability-pf transmission through 
the barrier varies substantially.ll2l For relatively large and 
shallow dots, such as those that have been constructed-in, 
GaAs/AlGaAs heterostructures at low temperatures,trEl 
the energies Up and W tend to be much smaller than the 
Fermi energy Ep but much larger than both the single- 
particle level spacing S and the thermal energy k^T. As 
a result, such systems are jcharactcrized by the following 
hierarchy of energy scales:Ej 

kBT,S<$:Up<W <^Ef. (1) 

Because W is comparable in size to the intermediate en- 
ergy scale Up (which acts as a kind of "excitation en- 
ergy scale" with respect to independent-particle ener- 
gies), effects from the barrier's finite size can be signifi- 
cant and deserve investigation before they are confidently 
discarded. 

Our prior study of the weak-coupling limit where g <C 1 
found that the finiteness of the barrier leads to an up- 
ward correction to the universal f-versus-g curve for a 
delta-function barrier. In other words, for a given small 
value of g, the fractional peak splitting / is erJaanced 
relative to its value for a delta-function barrierJl3 This 
enhancement occurs because, in a channel containing a 
barrier with a finite energy height (as opposed to a delta- 
function barrier), electrons can tunnel from one dot to an- 
other through largely unreflected states that have single- 
particle energies greater thgitLthe energy that corresponds 
to the finite barrier's peak.tj 

In order to study the correction to the universal /- 
versus-(/ curve, in our prior work we made several basic 
assumptions about the system of two symmetric dots con- 
nected by a single orbital channel. First, we assumed that 
the electrons that enter the point contact between the 
dots are in the lowest energy eigenstate for motion per- 
pendicular to the channel and can therefore be described 
by a one-dimensional (ID) Schrodinger equation with an 
effective barrier potential V{x). We assumed that V{x) 
could be treated as parabolic, with a harmonic oscilla- 
tor frequency top associated with the inverted parabolic 
well. (In this case, the transmission energy scale W is 
given by fuvp /2Tr.) In addition, we assumed that, de- 
spite the nonzero width of the barrier, the essential na- 
ture of electron-electron interactions was still accurately 
represented by the standard Coulomb blockade model, in 
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which the interaction energy of a quantum dot is anal- 
ogous to that of a classical capacitor — i.e., proportional 
to the charge on the dot squared. Having made these as- 
sumptions, we found that, in the limit g — > 0, the leading 
behavior of the enhancement to / is roughly proportional 
to Nch{2TTUp/hujp)/\ \ng\. For the experimental systems 
with which we are concerned, where A^ch is 1 or 2 and 
where 2T:Up/hu;u^ ~ 1, the size of the enhancement is 
relatively small.11^ 

This paper extends study of such finite-barrier correc- 
tions to the strong-coupling limit where (1 — g) ^ 1. In 
doing so, we repeat the two assumptions described above. 
The assumption of a parabolic form fo r t J ae . imposed bar- 
rier potential is likely as good as before.ElJE3 However, our 
simplistic treatment of the electron-electron interactions 
is likely less firm, given that the point contact is now 
largely open and therefore much more likely to be occu- 
pied by electrons interacting in a way not reflected by the 
standard capacitive Coulomb blockade model. Nonethe- 
less, because the barrier region is still small in relation to 
the larger conducting basins of the two dots, and because 
the effects from a combination of the Coulomb blockade 
model and the interdot barrier's finite size are themselves 
sufficiently interesting and complex, in this initial stab 
at the problem of finite barriers we ignore the effects of 
interactions specific to the point contact, with the under- 
standing that separate efforts to understand those effects 
should follow. 

Thus, the basic parameters for our strong-coupling 
study are the same as those for the weak-coupling limit. 
The Nch degenerate spin channels are assumed to contain 
effectively identical barriers, and the interdot barrier can 
be modeled as fully parabolic so long as we are concerned 
only with energies within a restricted range about the 
barrier peak. We therefore describe the single interdot 
orbital channel as containing the potential: 



V{x) = Fo (1 - x^/2e) 
V{x) = 



for \x\<V2^ 

for |a;| > V2^ . (2) 



See Figure 1. As we will later see, our ultimate results 
are independent of the potential's details away from the 
barrier peak. Key aspects of the barrier potential are the 
"depth" (Ep — Vo) of its peak relative to the Fermi surface 
and the harmonic oscillator frequency ujp ^ \/Vo/m^'^ 
associated with the inverted parabolic well.E£l (Note that 
fkup/Ep fa \/2 {kp£_)~^ , which is assumed to be much 
less than 1.) 

We use the single-particle eigenfunctions that corre- 
spond to the above potential to calculate corrections to 
/. To do so, we incorporate these eigenfjipctions in a gen- 
eralization of a bosonization approachEj Ijpthat was used 
earlier in the delta- function barrier limit We first find 
the change in the two-dot system's ground state energy 
as the value of (1 — g) is increased from zero. This energy 
shift A(p) is a function of {1 — g) and another dimension- 
less quantity p, which represents a linear combination of 
the gate voltages applied to the two dots.B'0 When the 



total number of electrons on the two dots is even, the 
fractional peak splitting / is given by the following for- 
mula: 



/ 



A(1)-A(0) 



(3) 



By applying this fornmla, we find that, consistent with 
our earlier conjecture,E3 for both A^ch ~ 1 and iVch — 2 
the leading finite-barrier corrections decrease the value of 
/ observed for a given value of g when g is close to 1 . For 
2TTUp/hujp < 2 and {1 — g) <C 1, the magnitude of the 
decrease is proportional to {2TTUp/hujp)/\\n{l — g)\ for 

ln(l-:9)|-i} 



Nch = 1 and to {2TTUp/hLUp) ^/l 
for = 2. 




^ = 2"% 



FIG. 1. A schematic representation of the ID model for a 
parabolic barrier between two quantum dots. The half width 
of the barrier at half maximum is ^, and the potential energy 
at the barrier's peak is Vo, which is less than, but comparable 
in size to, the Fermi energy Ep- 

In Sec. II, we show how to use the bosonization ap- 
proach to derive an effective action that captures the 
Coulomb blockade behavior of such a system for an ar- 
bitrary number A^ch of interdot tunneling channels. In 
Sec. HI, we use the effective action of Sec. II to solve for 
the leading /-versus-g behavior for Nch = 1 and 5 ~ 1. 
In Sec. IV, we do the same for A"ch = 2. In Sec. V, we 
remark on the significance and limitations of our results. 



II. BOSONIZATION WITH A BARRIER OF 
FINITE WIDTH 

To calculate the corrections due to the interdot bar- 
rier's nonzero width, we start with a fermionic formu- 
lation of our system of two symmetric quantum dots 
connected by a single orbital channel containing an ar- 
bitrary number Nch of degenerate spin modes. Since, 
this fermionic problem is effectively one-dimensional ,E2I 
we can apply toJit a variation of the standard technique 
of bosonizationJJ'B Having "bosonized" the problem, we 
integrate out various bosonic degrees of freedom to ob- 
tain a low-energy effective action involving only a single 
set of scalar fields. 

The fermionic formulation of the pfioWeiHi of coupled 
quantum dots is quite straightforward.Qfli3ia In a basis of 
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non-interacting single-particle eigenstates, the fermionic 
Hamiltonian consists of a kinetic-energy part Hk and a 
charging-energy part He'- 



(4) 



where ^fc is the energy of the single-particle eigenstate; a 
is the channel index; C is the index for states incident on 
the barrier from the left (C = 1) and the right (C = —1), 
respectively; and kp + k is the magnitude of the wave 
vector associated with motion along the one-dimensional 
channel (kp is the Fermi wave vector and |fc| is assumed 
to be smaller than kp). The operator h measures the dif- 
ference between the number of electrons to the right of 
the barrier and the number of electrons to the left of the 
barrier. If one lets hi be the occupation number for the 
left dot and n2 be the occupation number for the right 
dot, h = (n2 — hi)/2. In order to determine the value of 
/ via Eq. |[ we consider the case where the total number 
of electrons on the two dots is even, so that fi is an op- 
erator with integer eigenvalues. (When the total number 
of electrons is odd, n has half-integer eigenvalues.) 

The operator h can be written in terms of fermionic 
position operators '^'^{x) and ipa-{x) as follows: 



(5) 



where 4'a-{x) consists of left-incident and right-incident 
parts, 



(6) 



and where 0(x) is a Heaviside step function centered on 
the dividing line between the "right" and "left" sides of 
the barrier at a; = 0. 

If one were bosonizing fermionic operators associated 
with a basis of plane- wave states, one could bosonize the 
above Hamiltonian by using the standard formula: 



■Ipacix) = 



1 



(7) 



where Oa-{x) and (paix) are scalar boson fields and a is 
a quantity analogous to a lattice spacing that is used 
to impose an ultraviolet cutoff at energies aperpximat- 
ing hvp/a on the bosonic degrees of freedornllxilS (a is 
ultimately meant to be taken to zero — a limit that can 
be reasonably taken when Ep ^ Up and the approach to 
modeling the double-dot system can be supposed to work 
down to length scales a, where hvp/a is much greater 
than the relevant excitation scale Up). The annihilation 
operators associated with the plane-wave states would 
then be given by the equation: 
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O-kaC — 



dx e"^^^'' e*^'^^''^^''"'"'^"*-^-" (8) 



In terms of this plane-wave bosonization scheme, the op- 
erate* h could be simply expressed in terms of the boson 



^EMO). 



(9) 



Unfortunately, our problem is not so simple. The ba- 
sis of single-particle states that appears in Eq. ^ is not 
a plane-wave basis. The eigenstates are those of a one- 
dimensional system containing a barrier of finite height 
and nonzero width, and they therefore include incident, 
reflected, and transmitted parts. 

Happily, the kinetic-energy part of the Hamiltonian 
Hk does not concern itself with such complications. So 
long as the single-particle level spacing can be treated 
as effectively zero, we can proceed with regard to 
just as if the single-particle states were plane waves. Ac- 
cordingly, if we focus on bosonizingTJ/f , we can bosonize 
the annihilation operators in Eq. in the same way as 
in Eq. |^. Under this bosonization scheme, the kinetic 
energy Hk has the following standard form: 

Hk = yT. f {[dM^)? + [d.O,i^)f} . (10) 



Our next task is to find a proper bosonized form for the 
position operators ^crC (^) that appear in Eq. |^. Recalling 
that we already have bosonized forms for the operators 
afeo-Ci we express the ipa-ci^) in bosonized form by us- 
ing the fact that the position operators consist of linear 
combinations of the operators a/co-^- In particular, in the 
continuum limit. 



'^Tci^) ^ dk Yfcf (a;) a^ac 



(11) 



where the functions Yk(^{x) are the single-particle eigen- 
functions of the non-interacting one-dimensional sys- 
temrj For the potential described by Eq. |^, we have de- 
rived these Yk({x) in prior work, using the fact that there 
are exact solutions (parabolic cylinder functions) to the 
Schrodingcr equation for a single particle in a parabolic 
potential. t3 Using the yfc^(x), we find that, to leading 
order in the perturbation due to the barrier, n equals 
fiQ -\~ dn, where 



^^^I^ + C.C.J a,^,,a,,.c 
Sh = ^ ^ Jdki Jdk2 i?(ei, €2) 



^i[D{e2)~D{ei)] 

k2 - ki + irj 



4ial«fc2a,-l + h-c. \ , (12) 
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where "c.c." and "h.c." stand, respectively, for the com- 
plex conjugate and hermitian conjugate of the preceding 
term, and the following identities hold: 



i?(ei,e2) 



Die) 



2VT- 



-27re2 



argr(l/2-ie)-e + eln|e| 



+ C, (13) 



with C being a constant independent of e and^^ = 
^^fiL^^"'' being the value of e at the Fermi energy.o 

In the above equations, the dimensionless variables 
^« ^ ^'/iJ^^^^ measure the energy distance from the fi- 
nite barrier's peak.t2l The scale for these dimensionless 
energy measures is the harmonic oscillator energy hiup, 
and in terms of these dimensionless measures 
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1 



1 



(14) 



The functions R{ei,e2) and D{ei) are associated with 
the single-particle transmission properties of the barrier. 
Each additive term in R{ei, £2) is a product of transmis- 
sion and reflection amplitudes, reflecting the fact that 
the leading contributions to 6h come from "overlaps" 
between the transmitted part of an original right-mover 
and the reflected part of an original left-mover and vice 
versa). Meanwhile, D{ei) represents a phase associated 
with the scattering of a single-particle wavefunction. Be- 
cause of the exponential growth (or damping) of i?(ei, £2), 
we should not be surprised to find that i?(ei, £2) plays a 
starring role in the nature of our leading-order result — 
whereas, to first approximation, contributions from D{ti) 
will prove to be negligible. 

Before we proceed further, we should make three re- 
marks about Eq. First, we should note that derivation 
of the above equation for 5n ignores contributions to ng 
from integrating over some interior portion of the barrier 
region in which approximating the single-particle wave- 
functions by plane waves, or at least by WKB approxima- 
tions to plane waves, is no longer valid. Hence — because 
the exact eigenfunctions in the barrier region (parabolic 
cylinder functions) do not "explode" in ampilitude as one 
goes further toward the center of the barriered — as in our 
derivation of [analogous identities in the limit of weakly 
coupled dotsjH] we expect additive corrections of rough 



order ^ to the Eq. 02 cofactors 



Consequently, 



k2 — ki±iri' 

to avoid substantial corrections to Eq. [l^, we must be 
able to restrict our attention to wave-vectors ki (mea- 
sured from the appropriate Fermi wave vector in one di- 
mension) such that 



|fc.k<l. 



(15) 



Because the exponentially variant function i?(£i,£2) 
strongly favors values of ki that correspond to Ci near 0, 
the ki of most interest for calculating corrections to the 



universal behavior are those for which ki ~ —dpcf/vf- 
The above restriction on ki can thus be seen to mean that 
we need £f ^ \/'^Ep/Vo. For systems (such as those this 
paper seeks to model) in which Ep/V{) ~ 1, the restric- 
tion then reduces (conservatively speaking) to £i? ^ 1 or, 
equivalently, {\ — g)^ 10^"^. Given that|-|the relevant ex- 
periments with quantum dot systemaiTEl have generally 
not resolved g to increments of lO"'^, this constraint is 
not a serious one. 

As a second point regarding Eq. [ij, it should be noted 
that, for £j? < 1, this intermediate result (and therefore 
those that follow from it) is insensitive to details of the 
barrier away from the barrier peak (e.g., the slope discon- 
tinuities that occur in our barrier potential at x = ±-^2 ^) 
so long as ^ is large compared to the inverse Fermi wave 
vector k'p^ (with the Fermi wavelength itself being of the 
order of 2-Kh/ y/2mVo, the wavelength associated with the 
peak barrier energy, Vq). The condition ^ kp^ is typi- 
cally true of the kinds of ajiantum dots that have moti- 
vated these investigations Bi3'E£l 

A third remark about Eq. |l2| is that in deriving Sh, 
we have ignored a correction involving terms of the form 
(cos[iI'(£2) — D{€i)] — l)/{k2 — ki± IT]). We ignored sim- 
ilar "phase-based" terms in Ref. 10. The basic reason 
that we ignore these terms here derives from the fact 
that, because the characteristic excitation energy Up is 
less than hujp, the variation of the numerator of these 
terms over any region of particular interest should be 
very small. As a result, in the regions that we gener- 
ally find to be most important (those regions near the 
simple poles where /c2 = fci ± irj and the magnitudes 
of other contributing terms tend toward infinity), these 
phase-based terms go to zero. Thus, we have good reason 
to expect that the contribution from these phase-based 
terms is small compared to the contribution from the 
terms we keep. Indeed, if we follow the calculation of 
the lowest-order contribution from these terms further, 
we would find that, by reasoning parallel to that used in 
Sec. Ill to perform various wave vector integrals, these 
terms produce a contour integral with no contribution 
from residues of the simple poles at ^2 = fci ± it] and 
only negligible, higher-order contributions from integra- 
tion around other singularities and branch cuts. 

Returning to Eq. |l^ and the process of determining 
the fractional peak splitting, we observe that if the op- 
erators OfccrC were annihilation operators for particles in 
plane- wave states (e.g., if there were no interdot barrier), 
n would simply equal uq. In other words, no bears the 
same relation to the bosonic fields associated with the 
barrier-state operators akaC that n bears to the bosonic 
fields associated with the similarly indexed plane-wave 
operators of Eq. |[ From Eq. 0, it therefore follows that 



(16) 



We do not have a similarly simple formula for dn. Sub- 
stituting for the operators UkaC^ we find that 



5 



(2 



-^g— ^ Jdxi Jdx2 Jdki Jdk2 R{ei,e2) 



t,i[-D(£2) — -D(ei)]gi(fclKl+fc2 2;2) 

k2 - ki + irj 



h.c. 



(17) 



Having now obtained bosonized expressions for the 
separate terms fiQ and 5n, we express the Hamiltonian 
as a whole in bosonized form: 



H = H 



K 



H, 



(0) 



H, 



(1) 



H, 



(2) 



(18) 



Hk is given by Eq. and the remaining terms are as 
foUows: 



H, 



(0) _ 



H. 



(1) 



H. 



(2) 



Up Sh 



Up {Snf 



h.c. 



(19) 



(2) 

The term can be dropped from the leading-order 
calculation because, to second-order in dn (which due 
to i?(ei,e2) is itself roughly proportional to y/1 — g at 
the Fermi surface), Hq contributes nothing to the p- 
dependence of the ground-state energy. Recall from Eq. ^ 
that it is the p-dependence of the ground-state energy 
that provides thje_basis for calculation of the fractional 
peak-splitting /.ErO 

To simplify calculation of this p-dependence, we shift 
the 0cr-fields by a term linear in p: 9„{x) —> 9„{x) + 
y/T:p/2Nch_- This transformation leaves us with a Hamil- 
tonian in which the p-dependence appears only in the 
perturbative factors of 5fi. The transformed (and trun- 
cated) Hamiltonian has the following form: 

-1 2 



H. 



(0) 



Ur, 



IT 

u„ 



where dn is now given by the identity: 

= ^2^)30; E y y^2;2 Jdki J dk2 R{ei,e2) 



(20) 



„i[D(e2)--D(£i)]gi(fciXi + fe2a:2) 

k2 ~ ki + irj 

-i^/Tfle„{xi)+<t>„(xi)] 



X e 



X e 



-i-/TT[e^{x2) — 4>i7i^2)] 



h.c. 



(21) 



With Eq. we have reached the end of the road with 
regard to the Hamiltonian approach. To progress further, 
we shift to an action-based, path-integral formulation. 
For ease of notation, we will henceforth assume that we 
have chosen units in which 
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(22) 



and drop h from subsequent intermediate calculations 
(although we resuscitate it in stating our final results). 
Having made this choice of units, we proceed by inte- 
grating out various degrees of freedom: first, the 0-field 
degrees of freedom, and second, the 61- field degrees of free- 
dom away from x — 0. The result is an effective action 
in terms of fields 9^'' (t) that are equivalent to the 6'-field 
degrees of freedom at a; = 0. This effective action, depen- 
dent only on the a; = degrees of freedom for the original 
0-fields, is analogous, though — significantly — not equiva- 
lent, to that previously obtained for a system containing 
a delta- function barrier (see Refs. 1 and 2). 

We do no more than outline the process of integrat- 
ing out the 0-fields and 6'-fields because the methodol- 
ogy for doing so is fairly straightforward. With regard 
to integrating out the 0-fields, a key point is that, in 
the path-integral approach, a Lagrangian term linear in 
dx4>{x, t) dr9{x, t) appears in the action (just as a pi- 
term appears in going from the Hamiltonian for a single 
particle to the corresponding Lagrangian). Having ac- 
counted for this term, we eliminate the i/i-fields by per- 
forming a standard Gaussian integration. The result is 
an effective action entirely in terms of 0-fields. 

Having obtained this effective action, we introduce the 
field (t) and the auxiliary fields Ag- (r) . We use the A- 
fields to enforce the identity 6*^(0, t) ~ 9^\t) by adding 
to the action the following term: 



Sx 



E 



dT ^A,(r) [0^°)(r)-0,(O,T) 



(23) 



Having added 5a to the action, we can substitute 9'^'' (r) 
for 9^(0, t) in H^^ (recall Eq. |o|). 

Our next steps are to integrate out the 9c,{x,t)- 
fields and then the Acr(T)-fields. Once again, only 
straightforward — albeit tedious — Gaussian integrations 
are required. To leading order in the perturbative fac- 
tor Sh, the result is an effective action dependent only 
on the scalar fields 9^^\t) and their Fourier transforms 



(u;): 



E/S Ho'^\^)~ef\--) 



EE /£(^)^1?H^1"2H-) 



— ^ ^ jdr j dxi J dx2 J dki J dk2 
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2C/„ 



i?(ei,e2) A(xi,X2) 



X < e 



g-i[-D(e2)--D(ei)]g-i(feixi + fe2a;2) 



k2 — ki — ir] 



+ m.t.} , 



(24) 



where the "mirror term" (m.t.) can be obtained by com- 
plex conjugating both the factors that precede 9'§'^{t) 
and the exponential prefactor i that precedes the integral 
over w, and where the functions A{xi,X2), /i2(a;i, 2:2, cj), 
hi{xi,X2,ijj), and hQ{xi,X2,u}) are given by the following 
formulae: 



A{xi,X2) 



1 



(1^1 

X 



\X2\ + a 



,1/2 



(2|a;i|+a)i/4 (2|a;2| + a)i/4 



a{2\xi 



)i/2(2 



X2\+a 



,1/2 



1/2 



(|a;i| + a)(|a;2| + a) 
(|xi|+a)(|x2|+a)T^*^"(^^)^«"(^^)/' 



a(|a;i| + \x2\ + a) 

h2ixi,X2,Uj) = hi{xi,X2,Uj) - ho{xi,X2,U!) 

. , , sgn(x2) [1 - e-l""^l/-^] 
«i (a;i , a;2 , w) = sgn(tj) „ ^ ^ 

- {xi ^ X2} 

^~\u)Xi\/vf ^_ ^-\uJX2\/vf 

ho{xi,X2,Uj) = 



2^/^ 



(25) 



where "{xi <-^- X2}" is the same as the preceding term, 
except with the roles of xi and X2 reversed. 

The effective action of Eq. |4| is the end result for this 
section. As promised, this action depends only on a sin- 
gle set of scalar fields, which are equivalent to the original 
6'-fields at X = 0. Sections III and IV use this action to 
solve for the fractional peak splitting / in systems with 
one and two interdot channels, respectively. In both sec- 
tions, the general approach is to treat S'^"-' = + Sg 
as the unperturbed action and to solve perturbatively in 



III. FINITE-BARRIER RESULT FOR THE 
ONE-CHANNEL PROBLEM 

Here we consider the single- channel case (A^ch — 1), in 
which electrons of only one spin have enough energy to 
penetrate into the point contact. Starting with the un- 
perturbed action sf'^ — + , it is not hard to find 
that 



{9<-°\uj,)9'^"\lu2)) = 



^1 



C/p/TT 



d{ui+u;2). (26) 



We use this identity to calculate the first-order correc- 
tion to the ground state energy: A(p) = (S'g^'')//?. In 
particular, we find the following: 

A(p) = ^^^^^ jdxi jdx2 jdki jdk2 R{ei,e2) 

gC2{xi,X2,2Up)/'i 

{\xi\ + \X2\+ a)i/2 (2|a;i| + a)i/4 (2|:e2| + ay/^ 



X e 



-(1I-/2) / dui lhoixi,X2,uj)f 



X / dujho{xi,X2,uj) 



■]lu\/vf 



U)\ + Up/ IT 



-i[D(c2)~-D(iii)] ^-i(kiXi+k2X2) 

fc2 - fci - i-q 



+ c.c. 



(27) 



where the function C2(a;i, 2:2, 2J7p) is defined below in 

Before we apply standard techniques of complex anal- 
ysis to reduce the integral over k2 to the contributions 
from the simple poles at k2 — ki ± iij, we need to show 
that we can neglect the effects of the singularities in the 
expression for 5^^^ that result from the factors i?(ei,e2) 
and D{e). It is not too hard to discount the singular be- 
havior at £2 = that results from the functional form of 
D{e) (see Eq. O). We can rewrite eln |e| as follows: 



eln |e| = elne — i7reO(— e) . 



(28) 



where 0(— e) is a sort of variant of the Heaviside step 
function that, in the appropriate limits, is 1 for e real 
and negative and is for e real and positive. To allow for 
complex values of e, we can define 0(— e) by the formula 



6(-e) = lim^^o Jri^) dz 



r(€) 7T(z'i+v'i) 



, where the path r(e) 



in the complex plane starts at z = i? (with R being a 
large and positive real number, ultimately taken to 00) 
and then proceeds to the value e by the shortest path 
that avoids the singularities at z — ±iv and their accom- 
panying branch cuts. In the limit i? — s- cx), this definition 
yields the desired behavior for real values of e. On the 
other hand, it does leave Q{—e) with two branch points 
(at e — ±w), and to these we must add the branch point 
of elne at e = 0. 

Fortunately, it is not too hard to see that integrat- 
ing around these singularities and associated branch cuts 
leads to contributions to the final result that are higher- 
order in the expansion parameter 2i:Up/hiijp than the 
contributions from the simple poles at k2 — ki ± irj. 
Rough evaluation of the higher-order contributions indi- 
cates that they are indeed negligible in the regime where 
27TUp/hujp ^ 2, but are probably not negligible for values 
of 2TTUp/hujp approaching 10. 

Having discounted higher-order contributions from sin- 
gularities other than the simple poles at ^2 = fci ± ir], we 
can perform the k2 integral with ease. The integral over 
ki can then be done after the resulting factor i?(ei, ei) is 

expanded in powers of e^^"^^ (and after it is noted that 
g-i[_D(£i)-D(£i)] ^ g^j^j i-j^g j.gg^i^ jg ^jjg following: 
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2U 

^(P) = -JT^ cos{np)Ji{eF,Up,LOp) , (29) 

where, once again, e^^^ is the "dimensionless Fermi en- 
ergy" relative to the barrier peak, ujp \s the harmonic 
osciUator energy of the parabohc barrier (recaU the dis- 
cussion of Sec. I), and 

/OO POO 
dxi I dX2 [Q{xi) + e(-X2)] 
-OO J —OQ 

xe"^"^("i+"^)/"^Fi(a;i,a;2,C/p,c^p) . (30) 

As before, 8(2:) is the Heaviside step function and 
vp is the Fermi velocity. Furthermore, the function 
Fi(xi,X2,Up,LUp) has the formula: 



Fi{xi,X2,Up,LL!p) 



J- 



Up 



\2^/ttvp 



n=0 



2{2n + l)iT 



{2n^lYn-^ + {xi+X2YujI/vI 

^C2{xt.X2,Up)/i 

{\x,\ + \x2\+aY/^ 

1 

(2H+^FM2N+^)^ 

— {n/2) J duj [ho{xi,X2,ui)]'^ - 



^-a\uj\/vp 

X I duha{xi,X2,Lu) ———— , (31) 
\Lu\ + Up/n 



where 



C2ixi,X2,Up) = Ci{xi,X2,Up) + Ci{x2,Xi,Up) 



Cl{xi,X2, Up) 



dz ■ 



~{aUp/-nVF)z 







I _ q-{Up\xi\/t^vf)z 

- sgn(a;i) sgn(a;2) (l - e-^^-'l^^l/""^)^ 

X ^1 _ g-{Up\xi\/TTVF)z^ 

^^U,\x,\/^VF _ 1^ (1 _ sgn(2;i) sgn(a;2)) 

^-{a+\xi\)(Up/T^VF)z 

X dz 

Ji z 



X I dz 



■00 ^-{a+2\xi\)(Up/izvF)z 



■sgn(xi)sgn(a;2) (e'^'-d^^l+I^^D/'^^'^- - l) 



X I dz 



00 ^-{a+\xi\ + \x2\){Up/TrVF)z 



(32) 



Although the equation for A(p) may not be entirely 
transparent, it is not hard to confirm, by taking the 



relation ^/T^ g 
Ji{eF,Up,u!p) - 



limit cjp — > 00, that the result for / that follows from 
Eqs. ^ and |2| agrees lyrth that previously derived for a 
delta-function barrier flia T he key is to recognize that the 
, e^'^'^-f /Vl + e-'^-^^p means that that 
2e^-^l — g as 27rJ7p/a;p 0, where 7 
is the Euler-Mascheroni constant (7 ~ 0.577). 

The behavior for 2TTUp/hujp ^ can be found through 
numerical integration. The results for various values of 
2T:Up/hLup are displayed in Figure 2. As prjedicted in our 
study of the weak-coupling limit (g — > 0))13 for a given 
value of g in the vicinity of 1, the fractional peak splitting 
/ is reduced relative to that for a delta-function barrier. 
This strong-coupling depression of the /-versus-y curve 
becomes greater for larger values of 27rUp/hujp. 



1.00 



0.90 



0.80 
0.98 




0.985 



0.99 



0.995 



1.00 



FIG. 2. Plots of the leading (1 - 5) -> behavior of /, 
the fractional peak splitting, as a function of g, the dimen- 
sionless interdot channel conductance, for a one-channel con- 
nection between two quantum dots (A'^ch = !)■ Each curve 
corresponds to a different value of the quantity 2TiUp/luup 
(see legend to the left of the curves). The solid line is the 
result for an interdot barrier that has effectively zero width 
(2TiUp/hLjp = 0). The dashed and dot-dashed curves show 
the leading /-versus-g dependence for finite barriers with 
2-kU pjhLop taking on values from 0.5 to 32. The curves are 
expected to be substantially accurate at least for {\ — g) 
greater than or equal to approximately 10~^ and for values of 
2TiUp/hL0p not much larger than 2. 

For 2'!TUp/hu!p ~ 1, however, the downward correction 
is quite small. Since earlier approximations of the barrier 
as a delta function were designed to predict the be 
of experimental systems in which 2'KUp/hijjp ~ 1 
our calculation confirms that those prior results were sub- 
stantially correct. Figure 2 gives a sense of what the lead- 
ing 27rL'p/ftti;p-dependent corrections look like for small 
and large values of 2iTUp/hu!p, but, because of the ap- 
proximations we have made, the results are expected to 
be rehably accurate only for 2TTUp/hujp < 2. 

One might wonder what analytic form the leading 
27rC/p/?ia;p-dependent corrections take. In the weak- 
coupling limit {g —>■ 0), the finiteness of the barrier was 
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shown to result in leading-order correctioris-proportional 
to {2TTUp/niop)/\\a.g\ for 2T:Up/hwp < 2MCorrections 
of this form arose from the fact that transitions into vir- 
tual states with energies near or above that of the barrier 
peak allow for relatively free interdot movement (with 
transmission probabilities on the order of 1). 

By reasoning analogous (but in a sense opposite) to 
that used in studying the weakly coupled system, we ex- 
pect that, to leading order, the finite-barrier corrections 
in the strong-coupling limit, where {1 — g) <C 1, are dom- 
inated by the enhanced backscattering of states near or 
below the barrier peak, where reflection probabilities are 
on the order of 1, rather than on the order of ^/l — g as 
is true near the Fermi surface. To estimate these cor- 
rections from strongly backscattered states, it is best to 
return to the form of A(p) that appeared in Eq. |27 . 

Examination of the x-dependent factors in Eq. g7's in- 
tegrand reveals that the primary contributions for the 
integral come when the variables Xi have values of rough 
magnitude i:hvF /2Up. In this range of x^-values, the ki- 
independent factors of the integrand (i.e., all the factors 
but i?(ei, €2) and those enclosed in the final set of large 
brackets) vary relatively slowly, only changing substan- 
tially when one or both of the Xi values changes by an 
addend of order Trhvp/^Up. Consequently, if the wave 
vectors ki have magnitudes greater than 2Up/'KhvF, the 
fc^-independent factors can, to leading order, be treated 
as constants for Xi ~ ■KhvpJIU p^ and the Xi-dependence 
of the integrands in Eq. can be assumed to be domi- 
nated by the factor e.^^i^^^^+^i^'^'^ , 

The question then is whether the hi values that con- 
tribute most substantially to the 27rC/p/?icL;p-dependence 
are large enough (in magnitude) for this assumption to be 
justified. As we have reasoned above, the leading contri- 
butions to 27rC/p/fiwp-dependence are expected to come 
from ki such that energy of the corresponding eigenstate 
is near or below the barrier peak. In other words, the 
most significant wave vectors are expected to be ones 
such that Ep + hvpki < Vq. Consequently, we expect 
these ki to be negative and to have magnitudes approxi- 
mating fco = {ujp/vp)ep, where (given Eq. |lj) 



(33) 



As noted in Sec. II, the anticipated importance of neg- 
ative wave vectors of approximate magnitude fco (which 
yield values for equal or near to 0) means that we 
must have {1 — g) > 10~^ for our general approach 
to be valid. Our desire to gain a more compact ana- 
lytic approximation to the leading 27rf7p//ki;p-dependent 
behavior now leads us to irnpose another requirement: 
fco ^ 2Up/nhvp. From Eq. ^3|, this requirement means 
that we need (I-5) < f,-{-2./^){-2.T^Up/hujp) _ Combining this 
upper bound on {1—g) with our previous lower bound, we 
find that our general approach and the new assumption 
that we propose to make (that the e'''*^'^^'^^"'"''^^^^-factors 
dominate the a^j-dependence of the integrands) are both 



valid only if IQ-^ < (1-g) < e-^-i/^){-i^u,/hup) ^ rpj^g 
assumption therefore allows a useful analytic approxima- 
tion (i.e., an approximation with a significant range of g 
values in which it is valid) when 2-KUp/fiwp is not much 
greater than 2. For 2'KUp/fvjjp values slightly over 10, on 
the other hand, the assumption's utility "breaks down" 
entirely, as the room between the upper and lower limits 
on (1— 5) vanishes. Nonetheless, because we have already 
restricted our attention to 2nUp/fkvp < 2, this "break- 
down" is of no concern; for the values of2TTUp/hujp of real 
interest, our new assumption for the purpose of a com- 
pact analytic approximation imposes no more restraint 
than our initial assumption that (1 — g) <^ 1. 

Having determined that we can assume e='='('^i^i+'^2a;2)_ 
dominance, we can now proceed with estimating the 
corrections due to wave vectors in the vicinity of fco. 
Approximation of the integrations based on the above 
assumptions leads to the conclusion that, to first ap- 
proximation, the correction to A(p) is proportional to 
{2nU p/hwp) /\ ln(l — g)\. Thus, the leading finite-barrier 
correction to the fractional peak splitting / may be 
roughly written as —ci{2T:Up/hjjp)/\ In (1 — g)\, where ci 
is a small positive number. The form of this correction 
is parallel to that found in the weak-coupling limit; in 
both cases, the leading 27r/7p/fia;p-dependent terms are 
proportional to i2iTUp/tiujp)/\\-!i{^^)\. The results dis- 
played in Figure 2 bear out our approximation to the 
analytic behavior of the leading strong-coupling correc- 
tions. The graphed corrections for < {2'KUp/hwp) < 2 
agree with the predicted analytic behavior, producing er- 
rors of less than about 10% for the behavior as a function 
of (1 — g) and of less than about 20% for the behavior 
as a function of 2nUp/fvjjp. The graphed values suggest 
that the value of ci is slightly less than 0.05. 



IV. FINITE-BARRIER RESULT FOR THE 
TWO-CHANNEL PROBLEM 



The downward shift of the /-versus-g curve for one- 
channel systems (iVch = 1) is paralleled by a similar 
downward shift of the /-versus-g curve for two-channel 
systems (A'ch = 2). However, the route to the two- 
channel result is less straightforward than the already 
somewhat tortuous route taken for the single-channel 
problem, largely due to the fact that we now have to 
deal with two 6^^-?ie\As. in Eq. rather than just one. 
The first step in dealing with this set of fields is to trans- 
form them to a more manageable pair — the "charge" and 
"spin" fields 9c and 6s, respectively — which are linear 
combinations of the 61 and O2 fields: 6c 



and 6s = 6f^ — 62'' ■ The bosonic action S then consists 
of the sum of separate charge and spin contributions to 
the unperturbed action So, and a perturbative term Si, 
which depends on both the charge and the spin fields. In 

(C) (S) 

particular, Sq = Sq + Sq and S = Sq + Si, where 



a(0) 



9(0) 



9 



duj lull ~ 



Si - 



2U„ 



dr / dxi / dx2 / dfci / (ife2 



X i?(ei, £2)^(2:1, 2:2) 

„-i[D(e2)--D(ei)] -i(feixi+fe2 2;2) 



k2- ki -ir] 
+ m.t.} 

X |g-(V2)/<i'^'i2(a;i,X2,uj)e""^^ es(i^) ^ ni.t.j , (34) 

where the "mirror terms" (m.t.) can be obtained in the 
same way as for Eq. ^ 

After integrating out the high-energy charge fields, one 
is left with an effective spin action 
where S^^^ is unchanged from Eq. ^ and where 



(S) 



dr I dxi I dx2 / dki / dk2 



(2^)30F Jo 
X R{ei,e2) A{xuX2) eC=^(^i.=^=^'2C/p)/8 

a(2|a;iH-a)i/2 {2\x2\ + a^/^^''^^ 



{\xi\ + a) {\x2\ + a) 

N2|+«)(Nl|+«)^^''"^"'^''"'"'^/' 



a {\xi\ + \x2\+ a) 

{n/4) Jduj[haixi,X2,uj)] ^^^^2Up/ 
^~-a\Lu\/vp 



X e 

y. I duo ho{xi,X2,Lo) 



\uj\ + 2Up/n 



-i[D(e2)-D{ei)] -i{kiXi+k2X2) 

ie^pi^- + c.c. 

k2- ki - ir] 



(35) 



where, as was done in the case of the iVch = 1 prob- 
lem, only the lowest-order perturbative term has been 
retained in the effective action. 

To exploit the even-odd combinations that occur in 
the integrals over w, we rewrite the action in terms of 
the following "even" and "odd" fields: 



5(o) _ 



(s) 

The nonperturbative term Sq then takes the form 
1 



(36) 



qiS) _ 



4 ./o 2TT 



du! 



(37) 



After integrating out the 6'g^^-fields, we have a new ef- 
fective action in terms of the Og -fields. 

The "even" fields are not to be left alone, however. 
In anticipation of "refermionization," we insert into the 
effective action a set of dummy "odd" fields — free fields 
that are entirely decoupled from the "even" fields. The 
nonperturbative, kinetic energy terms of these new 9^^^- 

fields and 9g -fields add to give a total kinetic energy of 

the same form as that in Eq. This total is in turn 

rewritten in a form equivalent to that of Sq^'^ in Eq. |3j. 

Our newest, and final, effective action then consists of the 
(s) 

sum of this Sq -term and the leading perturbative term 
Sl, which was derived from integrating out the original 
"odd" fields: 



S{ 



(27r)3 



/3 



dr / dxi / dx2 / dki / dk2 



R{ei,e2)e 



C2{xi,X2,2Up)/8 



{\xi\ + \X2\ + a)i/2 (2|xi| + «)i/4 (21x21 + a)^^ 

-(7t/4) / duj [hoixi,X2,uj)f 



l/v 



X I duj ho{xi,X2, oj) 



^-a\uj\/vF 



g-i[D(e2)--D(ei)] ^-i(kiXi+k2X2) 



k2 — ki — ir] 

J duj ha(xi,X2,uj) e^'"^ Ss(aj) 



+ C.C. 



m.t. 



(38) 



As in prior work,Q'Ll we "refermionize" the bosonic ac- 
tion S = S*^^' + Si by identifying ^/>/(0, r) as a fermionic 
annihilation operator at the origin of a semi-infinite sys- 
tem, where 

g-(i/2) J du ho(xi,X2,(jj) e~"^^ 8si<^) 

V'/(0,r) = , (39) 



y^2TT a{xi,X2) 



where 



aixi,X2) = i2\xi\+a)^/\2\x2\+a)^/\\xi\ + \x2\+ay/\ 

(40) 

This refermionization formula may look peculiar because 
of its use of a normalization factor and an ultraviolet 
cutoff (embedded in ho{xi, X2,uj)) that depend on posi- 
tion variables, rather than mere constants (contrast the 
standard bosonic representation of fermionic position op- 
erators in Eq. ^ . These position-dependent factors work, 
in tandem with the scalar prefactor (i/2) for the expo- 
nentiated integral, to ensure the correct anticommuta- 
tion relations for -0/(0, r), and their use is expected to 
be unobjectionable so long as for all Xi of real interest 
the resulting ultraviolet cutoff is high enough to capture 
the behavior with which we are concerned. 

As in the single-channel problem, the Xi of real interest 
satisfy \xi\ ~ TThvF/2Up (as can be seen by doing the inte- 
grals over ki and obtaining a result with factors analogous 
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to those of Fi(xi, X2, Up, ujp) of Eq. ^l|). The ultraviolet 
cutoffs in Eq. ^ therefore characteristically correspond 
to energies approximating 2Up/iT — which is much greater 
than ksT, the characteristic energy for the unperturbed 
spin degrees of freedom, and (as we can confirm once we 
rediscover the leading order behavior produced by S{), 
much greater than the energy of states characteristically 
brought into play by the perturbation S^. Indeed, we 
have already incorporated an assumption that ultravi- 
olet cutoffs approximating 2Up/TT are valid because, in 
obtaining our effective action, we have only kept the lead- 
ing order terms from integrating out the charge degrees 
of freedom — an approximation only expected to be good 
if the spin-field states of cancern have excitation energies 
significantly less than UpB 

Returning to Eq. we use the identity 

V2^J„^dT [V/(0,r) + ^}(0,t)] = Jdk ifl + fk) to find 

that our refermionization scheme produces a fermionic 
Hamiltonian H — Hq + Hi that consists of the following 
parts: 



Ha = Jdk flfk 

Hi = Z{eF,Up,p) jdk ifl+h), 



(41) 



where 

Z{eF,Up,up,p) 



Uo 



X 



(27r)3V^ 
i?(ei,e2)e'='=^(=^i'^^'2f^'')/« 



dxi / dx2 / dki / dk2 



i\xi\ + \X2\ + ay/^ i2\xi\ + a)i/8 (2\x2\ + a)i/8 

^-a\u)\/vF 



X / dio ho{xi, X2,i-^ 

p/2 



X 



w| -I- 2Up/n 



I e 



k2- ki - ii] 



(42) 



From prior work on the effects of a delta-function 
barrier, we know how to solve for the fractioaai peak 
splitting that such a Hamiltonian produces.LTEl The 
only difference from the delta-function barrier Hamil- 
tonian that we encountered before is that a com- 
plicated prefactor Z{eF,Up,Ljjp, p) replaces the sim- 
pler delta- function barrier prefactor Zoo{^F,Up, p) — 

cos(7rp/2) -y/1 — g yj2e^hvFU pji^'^ , d^M^'-^ ^ Euler- 
Mascheroni constant (7 ~ 0.577) J3~0 Thus, to find the 
leading behavior of the fractional peak splitting as a func- 
tion of (1 — g) and 27:11 p/hujp, we can simply substitute 
Z{eF,Up,LUp, p) for ZooiiF,Up, p) in the results previ- 
ously obtained for a delta- function barrier. 

As for one-channel systems, we ultimately resort to nu- 
merical integration to solve for the fractional peak split- 
ting when the barrier has a nonzero width. Needless 
to say, it is gratifying that such numerical calculation 
confirms that, at least through five significant digits, the 



prefactor Z{eF,Up,ujp, p) converges to the delta- function 
quantity Zaoi^p, Up, p) in the limit 2TTUp/hiOp — > 0. Con- 
sequently, in the limit of a narrow barrier, we recover the 
same result as for a delta-function barrienJEl — a good 
confirmation both of the robustness of our earlier results 
and of our success in wending through the complications 
created by the initial assumption of a nonzero-width bar- 



The approach to performing the quadruple integral of 
Eq. ^ is very similar to that used to perform the anal- 
ogous quadruple integral for Nch = 1 and therefore will 
not be described in depth. We first integrate over ki and 
k2 by techniques of complex analysis similar to those used 
for iVch = 1, once again focusing on the simple poles at 
k2 — ki ± ii], with the understanding that this imposes 
a constraint that 2'iTUp/fiiijp ^ 2. For various values of 
2nUp/hujp, we then do the remaining integrals over xi 
and X2 numerically. The results are shown in Figure 3. 
As for one-channel systems, we see that as 2'iTUp/fiujp in- 
creases, the strong-coupling end of the /-versus-^ curve 
shifts downward from the zero-width result. Once again, 
the corrections for the experimentally realized values of 

the 
was 



2'KUp/hhjp ~ 1 are small, a fact which confirms 
previous assumption of a delta-function barrie: 
substantially justified, at least so long as interaction ef- 
fects peculiar to the barrier region can be ignored. 



1.00 




0.985 



FIG. 3. Plots of the leading (1 — gr) — > behavior of /, 
the fractional peak splitting, as a function of g, the dimen- 
sionless interdot channel conductance, for a two-channel con- 
nection between two quantum dots (A^ch = 2). Each curve 
corresponds to a different value of the quantity 2TiUp/hwp 
(see legend to the left of the curves). The solid line is the 
result for an interdot barrier that has effectively zero width 
{2TiUp/hiup = 0). The dashed and dot-dashed curves show 
the leading /-versus-y dependence for finite barriers with 
2-kU pjTujp taking on values from 0.5 to 32. The curves are 
expected to be substantially accurate at least for (1 — g) 
greater than or equal to approximately 10~^ and for values of 
2-kU pjhjjjp not much larger than 2. 
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As for the single-channel system, one might inquire 
about the nature of the leading analytic behavior of these 
corrections. The answer is essentially parallel, both in 
reasoning and substance, to that found at the end of 
Sec. III. The constraints on (1 — g) are the same (i.e., 
10-3 < (1 - g) < e-(2/^)(2^t/p/^p))^ aj^j ^Yic leading 

finite-barrier corrections to Z{eF,Up, ujp, p) are expected 
to be roughly proportional to {iTiU p/hiop) /\h\{l — g)\ 
when 2TrUp/hu!p ^ 2. In other words, Z{ep,Up,ujp, p) 
can be expanded as follows: 

Z{eF,Up,UJp,p) ^ n^— 

cos{'Kp/2)^2e'^hvFUp/T:^ ^ 

27111 p/hiop 



C2 



ln(l-g)| 



(43) 



where C2 is a srnall positive number. Because, 
from prior workj3Ll the leading corrections to 
the fractional peak splitting are proportional to 
Z{ep,Up,ojp, p)'^ \n Z{eF,Up,ujp, p) aX p = 0, the lead- 
ing finite-barrier correction to the fractional peak split- 
ting is roughly proportional to {2i:Up/fvjjp) y/\ — g {1 — 
I In (1 — g)!"^}. The corrections shown in Figure 3 follow 
this predicted behavior quite well: for 2'KUp/hujp < 2, 
the analytic prediction captures the calculated correc- 
tions with error margins of about 10% or less for the 
(1 — g)-dependence and of about 20% or less for the 
27rJ7p/fiii;p-dependence. The numerical results suggest 
that C2 is slightly less than 0.02. 



V. CONCLUSION 

This paper shows that bosonization techniques for 
studying the behavior of one-dimensional systems need 
not be abandoned when finite-length barriers (or con- 
strictions) are introduced. The finite-length effects of 
those barriers can be captured to leading order by a 
standard perturbative approach, albeit one that requires 
complicated calculations and a fair degree of care. Thus, 
bosonization techniques can still be useful when nontriv- 
ial behavior in the barrier region is at issue, and the 
approach presented in this paper may help investigators 
to distinguish between single-particle and many-particle 
effects from a barrier's finite length. 

With regard to the more particular problem of the 
Coulomb blockade behavior of coupled quantum dots, our 
results can be summarized as follows. This paper shows 
that at least for one-channel or two-channel systems, the 
fractional peak splitting / of two strongly coupled dots 
decreases, for a given value of the interdot conductance, 
as the ratio 2TTUp/hujp is increased. For one-channel sys- 
tems (iVch = 1), the downward correction behaves as 
{2nU p/fvjjp) /\ h\{l — g)\ in the limit where 2TiUp/}Tiop < 2 
and (1 — g) <C 1. Thus, for iVch = 1, the fractional peak 
splitting has the following leading-order functional form: 



where ci,2 is gpHiewhat less than 0.05 and, as derived 
in prior workJ3~cl Ci^i — 8e^/7r^ (about 1.44). In two- 
channel systems {N^h = 2), and in the same limits, the 
downward correction due to the finite barrier width be- 
haves as {2ttU p / hhj p) ^1 ~ g {1 — | In (1 — g)|~^}, result- 
ing in leading-order behavior of the form 



/2-l-C2.i(l-5)|ln(l-.9)| 

- C2.2 {2^Up/fiwp) Vl -5 { 1 - 



ln(i-.g)| 



(45) 



'iWH 



/i = i-ci,iv/r^ 



Cl,2- 



(2TTUp/hujp) 

ln(l-5)l 



(44) 



where C2,2 is somewhat less than 0.04 and, as derived 
in prior work,LrH C2,i = IGe'^/Tr^ (or about 0.919). (It 
should be noted that prior results indicate that, in the 
A'ch = 2 equation, the next additive term independent of 
2nUp/hLUp will equal —€2,3(1 — g), where 02,3 — 0.425.0 
Because this term linear in (l — g) comes from higher- 
order terms in the effective action than those we have 
considered here (see, e.g., Eq. we have omitted it 
in our calculations and in the graphs of Figure 3 in or- 
der to have a truer comparison between the 2T:Up/hujp- 
dependent terms and those they directly correct. If need 
be, the linear term is easy enough to combine with the 
result expressed in Eq. [45| .) 

The above results combine with an earlier study of the 
weak-coupling regime {g <C 1)£3 to give a more complete 
understanding of the /-versus-g curve when one leaves 
the delta- function barrier limit 2TrUp/hujp — s- 0. In par- 
ticular, we come to the nontrivial conclusion that, fj 
experimentally realized values of 2'KUp/hujp 
the corrections to the results derived from modeling the 
barrier as a delta function are not fundamentally sub- 
stantial; thus, to this extent at least, earlier theoretical 
work using a delta- function potential was correct. 

In addition to confirming the essential nature of the 
/-versus-g curve for 2'KUp/hjjjp ^ 2, this paper has 
helped us gain a better picture of what happens to the 
/-versus-g curve for more general values of 2iTUp/hujp. 
The weak-coupling results suggested that, as 2TTUp/hujp 
is increased, the /-versus-g curve shifts upward for small 
values of g (i.e., for g <C 1) and becomes flatter for inter- 
mediate values of g. The strong-coupling results suggest 
(as conjectured) that, as the same ratio is increased, the 
/-versus-g curve shifts downward for large values of g 
(i.e., for (1 — g) <C 1) and becomes flatter for interme- 
diate values of g. Together, the two sets of results sug- 
gest that, in the limit of an extremely wide, "adiabatic" 
barrier {2i:Up/hjjp — > 00), the /-versus-g curve will be 
essentially flat for values of g in an intermediate region 
between and 1. For such an adiabatic barrier, the /- 
versus-g curve will sit at some essentially constant value 
of / for most of this interval and will turn sharply toward 
the limiting values of / = and / = 1 at the edges. 

For the moment, the intermediate value at which the 
"adiabatic" /-versus-g curve sits remains a mystery. The 
weak-coupling results displayed an antisymmetry around 
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(7=1/2 that, if true at higher orders, might aid in solv- 
ing for the fractioBa,! peak sphtting when g takes on in- 
termediate values.113 Unfortunately, the strong-coupling 
results do not exhibit such a simple symmetry. As Fig- 
ures 1 and 2 indicate, unlike the weak-coupling results, 
the strong-coupling results do not reveal a common point 
of intersection for the leading-order curves that corre- 
spond to different values of 2TTUp/hujp. The apparent 
lack of a common pivot will presumably make solution of 
the intermediate-^ problem more difhcult. 
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